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ABSTRACT 

The  purpose  of  this  study  is  to  investigate  compressibility  effects  on  the  turbulence  in  homo¬ 
geneous  shear  flow.  We  find  that  the  growth  of  the  turbulent  kinetic  energy  decreases  with 
increasing  Mach  number  -  a  phenomenon  which  is  similar  to  the  reduction  of  turbulent  ve¬ 
locity  intensities  observed  in  experiments  on  supersonic  free  shear  layers.  An  examination  of 
the  turbulent  energy  budget  shows  that  both  the  compressible  dissipation  and  the  pressure- 
dilatation  contribute  to  the  decrease  in  the  growth  of  kinetic  energy.  The  pressure-dilatation 
is  predominantly  negative  in  homogeneous  shear  flow,  in  contrast  to  its  predominantly  pos¬ 
itive  behavior  in  isotropic  turbulence.  The  different  signs  of  the  pressure-dilatation  are 
explained  by  theoretical  consideration  of  the  equations  for  the  pressure  variance  and  density 
variance.  We  obtained  previously  the  following  results  for  isotropic  turbulence;  first,  the 
normalized  compressible  dissipation  is  of  O(Mf),  and  second,  there  is  approximate  equipar- 
tition  between  the  kinetic  and  potential  energies  associated  with  the  fluctuating  compressible 
mode.  Both  these  results  have  now  been  substantiated  in  the  case  of  homogeneous  shear. 
The  dilatation  field  is  significantly  more  skewed  and  intermittent  than  the  vorticity  field. 
Strong  compressions  seem  to  be  more  likely  than  strong  expansions. 
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1  Introduction 


Homogeneous  shear  flow  refers  to  the  problem  of  spatially  homogeneous  turbulence  sustained 
by  a  parallel  mean  velocity  field  u  =  (5x2, 0,0)  with  a  constant  shear  rate  S.  Such  a  flow 
is  perhaps  the  simplest  idealization  of  turbulent  shear  flow  where  there  are  no  boundary 
effects,  and  where  the  given  mean  flow  is  unaffected  by  the  Reynolds  stresses.  Nevertheless, 
the  crucial  mechanisms  of  sustenance  of  turbulent  fluctuations  by  a  mean  velocity  gradient, 
and  the  energy  cascade  down  to  the  small  scales  of  motion  are  both  present  in  this  flow. 

The  homogeneous  shear  flow  problem  has  been  studied  experimentally  by  Champagne, 
Harris  and  Corrsin  (1970),  Harris,  Graham  and  Corrsin  (1977),  and  Tavoularis  and  Corrsin 
(1981)  among  others.  In  these  low-speed  experiments  the  statistical  properties  of  the  flow 
do  not  vary  spatially  in  the  transverse  ( x2 ,  X3)  plane  but  they  evolve  in  the  streamwise  Xj 
direction.  In  theory,  the  streamwise  inhomogeneity  can  be  removed  by  writing  the  equations 
of  motion  in  a  reference  frame  moving  with  the  mean  flow  u.  In  the  moving  reference  frame, 
the  one-point  moments  satisfy  pure  evolution  equations  in  time  clearly  illustrating  that  ho¬ 
mogeneous  turbulence  is  fundamentally  an  initial  value  problem.  Rogallo  (1981)  and  Rogers 
and  Moin  (1987)  have  investigated  the  incompressible  homogeneous  shear  problem  at  great 
depth  through  direct  numerical  simulations.  These  simulations,  albeit  at  low  turbulence 
Reynolds  numbers,  have  provided  turbulence  statistics  which  are  in  good  agreement  with 
experiments  performed  at  relatively  higher  Reynolds  numbers.  Furthermore,  since  the  sim¬ 
ulations  provide  global  instantaneous  fields,  the  turbulence  can  be  studied  in  much  greater 
detail  than  in  physical  experiments. 

Recently  there  has  been  a  spurt  of  activity  in  the  direct  numerical  simulation  (DNS)  of 
three-dimensional  compressible  turbulence.  Decaying  isotropic  turbulence  has  been  studied 
by  Passot  (1987),  Erlebacher  et  al.  (1990),  Sarkar  et  al.  (1989),  and  Lee,  Lele  and  Moin 
(1990).  The  simulations  of  Erlebacher  et  al.  (1990)  identified  different  transient  regimes 
including  a  regime  with  weak  shocks,  and  also  showed  that  a  velocity  field  which  is  initially 
solenoidal  can  develop  a  significant  dilatational  component  at  later  times.  Sarkar  et  al. 
(1989)  investigated  the  statistical  moments  associated  with  the  compressible  mode  in  their 
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simulations,  and  determined  a  quasi-equilibrium  in  these  statistics  for  moderate  turbulent 
Mach  numbers  which  was  then  used  to  model  various  dilatational  correlations.  Lee,  Lele 
and  Moin  (1990)  studied  eddy  shocklets  which  developed  in  their  simulations  when  the 
initial  turbulent  Mach  number  was  sufficiently  high  ( Mt  >  0.6).  Kida  and  Orszag  (1990) 
primarily  studied  power  spectra,  and  energy  transfer  mechanisms  between  the  solenoidal  and 
dilatational  components  of  the  velocity  in  their  simulations  of  forced  isotropic  turbulence. 

Physical  experiments  have  not  been  and  perhaps  cannot  be  performed  for  homogeneous 
shear  flows  at  flow  speeds  which  are  sufficiently  high  to  introduce  compressible  effects  on  the 
turbulence.  However,  direct  numerical  simulation  of  this  problem  could  provide  meaningful 
data,  especially  since  DNS  for  the  incompressible  problem  has  been  successful  in  giving 
realistic  flow  fields.  The  compressible  problem  was  considered  by  Feiereisen  et  al.  (1982) 
who  performed  relatively  low  resolution  643  simulations  and  concluded  that  compressibility 
effects  are  small.  Recently  Blaisdell  (1990)  has  aiso  considered  compressible  shear  flow. 
We  have  performed  both  963  and  1283  simulations  which  have  allowed  us  to  obtain  some 
interesting  new  results  regarding  the  influence  of  compressibility  on  the  turbulence.  In 
contrast  to  the  results  of  Feiereisen  et  al.,  our  simulations  which  start  with  incompressible 
initial  data  develop  significant  rms  levels  of  dilatational  velocity  and  density.  We  find  that 
the  growth  rate  of  the  kinetic  energy  decreases  with  increasing  Mach  number  as  well  as 
increasing  rms  density  fluctuations  and  show  that  the  compressible  dissipation  and  pressure- 
dilatation  contribute  to  this  effect.  Apart  from  rms  levels  of  the  fluctuating  variables,  we 
have  also  examined  their  probability  density  functions  (pdf)  and  higher  order  moments.  The 
dilatational  field  has  different  skewness  and  flatness  characteristics  relative  to  the  vorticity 
field.  In  particular,  the  dilatational  field  has  a  significant  negative  skewness  and  is  more 
intermittent  than  the  vorticity  field. 

2  Simulation  Method 

The  compressible  Navier-Stokes  equations  are  written  in  a  frame  of  reference  moving  with  the 
mean  flow  tq.  This  transformation,  which  was  introduced  by  Rogallo  (1981)  for  incompress- 
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ible  homogeneous  shear,  removes  the  explicit  dependence  on  u(x2)  in  the  exact  equations  for 
the  fluctuating  velocity,  thus  allowing  the  imposition  of  periodic  boundary  conditions  in  the 
x2  direction.  The  relation  between  x*  and  the  lab  frame  Xi  is 

Xj  —  X\  Six  2  ,  X2  —  X2  ,  X3  —  X3 

Here  S  .enotes  the  constant  shear  rate  u12.  In  the  transformed  frame  x*,  the  compressible 
Navier-Stokes  equations  take  the  following  form 

dtp  +  (pui').i  -  St(pu2'),  1  =  0  (1) 

dt{put')  +  {pu/ui)j  =  -pti  +  r'ijj  -  Spu2'6ti 

+  5t(pu2,ui,)il  +  StpASl2  -  Str'i 21  (2) 

dtp  +  u/P'j  +  7  pu'jj  =  S  tu2'ptl  +  7-S"  tpu2i 

+  (7  -  1  )K[Tjj  -  2 StTn  +  S2t2Tn]  +  (3) 

p  =  pRT  (4) 

where  $  =  rtJUij  is  the  dissipation  function,  ux'  the  fluctuating  velocity,  p  the  instanta¬ 
neous  density,  p  the  pressure,  T  the  temperature,  R  the  gas  constant,  ^.nd  k  the  thermal 
conductivity.  The  viscous  stress  is 

2 

r >j  =  +  •Uj.,)  -  -p{uktk)6tj 

where  p,  is  the  molecular  viscosity  which  is  taken  to  be  constant.  All  the  derivatives  in  the 
above  system  are  evaluated  with  respect  to  the  transformed  coordinates  x*. 

Since  Eqs.  (l)-(4)  do  not  have  any  explicit  dependence  on  the  spatial  coordinates  x*,  and 
because  the  homogeneous  shear  flow  problem,  by  definition,  does  not  have  any  boundary 
effects,  periodic  boundary  conditions  are  allowable  on  all  the  faces  of  the  computational  box. 
Of  course,  in  order  to  obtain  realistic  turbulence  fields  it  is  necessary  that  the  length  of  the 


3 


computational  domain  be  much  larger  than  the  integral  length  scale  of  the  turbulence.  Spec¬ 
tral  accuracy  is  obtained  by  using  a  Fourier  collocation  method  for  the  spatial  discretization 
of  the  governing  equations.  FFT’s  are  used  to  obtain  the  Fourier  representation  from  the 
data  in  physical  space  and  thereby  calculate  derivatives.  In  order  to  avoid  expensive  evalua¬ 
tions  of  convolutions  in  Fourier  space,  the  nonlinear  terms  are  directly  evaluated  in  physical 
space  as  products  of  derivatives.  It  has  been  shown  by  Canuto,  Hussaini,  Quarteroni  and 
Zang  (1988),  that  the  ensuing  aliasing  error  is  negligible  if  the  significant  spatial  scales  of  the 
computed  variable  are  resolved  on  the  grid.  A  third  order,  low  storage  Runge-Kutta  scheme 
is  used  for  advancing  the  solution  in  time. 

The  computational  mesh  becomes  progressively  skewed  because  the  mean  flow  velocity 
Uj  with  which  the  mesh  moves  has  a  variation  in  the  x2  direction.  In  order  to  control  the 
loss  of  accuracy  caused  by  excessive  skewness  of  the  mesh,  we  interpolate  the  soi  tion  on 
the  grid  at  St  =  0.5  shown  in  Fig.  2a,  onto  the  grid  sketched  in  Fig.  2b  which  is  skewed 
in  the  opposite  direction  and  corresponds  to  St  =  —0.5.  Since  the  solution  is  periodic,  the 
interpolation  is  straightforward  in  physical  space.  However,  the  remesh  procedure  introduces 
spectral  aliasing  errors  as  pointed  out  by  Rogallo  (1981).  In  the  dealiasing  procedure  used 
here,  the  higher  Fourier  modes  are  truncated  both  before  and  after  the  remesh. 

Initial  conditions  have  to  be  prescribed  for  ut',  p ,  p  and  T.  The  initial  velocity  field  is 
split  into  two  independent  components,  that  is,  ut'  =  u['  +  uf',  each  component  having  a 
zero  average.  The  solenoidal  velocity  field  u['  which  satisfies  V.u1'  =  0  is  chosen  to  be  a 
random  Gaussian  field  with  the  power  spectrum 

E(k)  =  k*  exp(-2 k2/k2J  (5) 

where  km  denotes  the  wave  number  corresponding  to  the  peak  of  the  power  spectrum.  The 
compressible  velocity  uf  which  satisfies  V  x  uc>  =  0  is  also  chosen  to  be  a  random  Gaussian 
field  satisfying  the  same  power  spectrum,  Eq.  (5).  The  power  spectra  of  the  two  velocity 
components  are  scaled  so  as  to  obtain  a  prescribed  uTTnt  =  J u'u',  and  a  prescribed  x  = 
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urm»/urm»  which  is  the  compressible  fraction  of  kinetic  energy.  The  pressure  p1'  associated 
with  the  incompressible  velocity  is  evaluated  from  the  Poisson  equation 


vy '  =  -2pSu''u  - 


(6) 


It  remains  to  specify  the  thermodynamic  variables.  The  mean  density  p  is  chosen  equal  to 
unity,  and  p  is  chosen  so  as  to  obtain  a  prescribed  Mach  number  u-n^l yj^pj'p  characterizing 
the  turbulence.  The  fluctuating  density  p'  and  compressible  pressure  pc>  are  chosen  as 
random  fields  with  the  power  spectrum  Eq.  (5)  and  prescribed  p^,  and  p^.  The  pressure 
then  becomes  p  —  p  +  p1  +  pc  ,  the  density  is  p  =  p  +  p',  and  the  temperature  T  is  obtained 
from  the  equation  of  state  p  =  pRT. 


3  Results 

We  have  performed  simulations  for  a  variety  of  initial  conditions  and  obtained  turbulent  fields 
with  Taylor  microscale  Reynolds  numbers  Re *  up  to  35  and  turbulent  Mach  numbers  Mt  up 
to  0.6.  Note  that  Re\  =  q\/v  where  q  =  yj u'tu't  and  A  =  q/ \J w'x lu-  ,  while  Mt  =  q/c  where  c 
is  the  mean  speed  of  sound.  The  computational  domain  is  a  cube  with  side  2ir.  The  results 
discussed  here  were  obtained  with  a  uniform  963  mesh  overlaying  the  computational  domain. 
The  initial  values  of  some  of  the  important  non-dimensional  parameters  associated  with  the 
DNS  cases  ranged  as  follows:  3.6  <  ( SK/e)0  <  7.2,  16  <  Re\, o  <  24,  0.2  <  Mfi0  <  0.4,  and 
0  <  x  <  0.15.  We  note  that  K  denotes  the  turbulent  kinetic  energy  and  e  the  turbulent 
dissipation  rate. 

In  previous  work  by  Erlebacher  et  al.  (1990)  on  isotropic  turbulence,  it  was  found  that  for 
a  given  initial  Re\,  the  initial  choice  of  the  non-dimensional  quantities  Mt ,  Pmu/p,  Tn^/T, 
Prm*/p  and  x  had  a  strong  and  lasting  influence  on  the  temporal  evolution  of  the  turbulence. 
The  different  choices  of  initial  conditions  that  lead  to  different  regimes  were  classified  by 
Erlebacher  et  al.  Even  though  the  long  time  statistics  in  isotropic  turbulence  are,  in  general, 
strongly  dependent  on  the  initial  conditions,  there  is  a  quantity  in  isotropic  turbulence  which 
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has  a  large  degree  of  universality  (in  the  absence  of  shocks).  This  quantity  is  the  partition 
factor  F  defined  by 

F  =  7  2M?x/p£  (?) 

which  was  shown  by  Sarkar  et  al.  (1989)  to  equal  unity  (to  lowest  order)  in  a  low  Mt 
asymptotic  analysis,  and  to  approach  unity  from  a  variety  of  initial  conditions  in  direct 
simulations  of  isotropic  turbulence  with  Mt  up  to  0.6.  The  physical  interpretation  of  F  — >  1 
is  that  there  is  a  tendency  towards  equipartition  between  the  kinetic  energy  and  potential 
energy  of  the  compressible  component  of  the  turbulence. 

In  the  case  of  initially  isotropic  turbulence  subjected  to  homogeneous  shear,  for  a  given 
initial  energy  spectrum,  the  solution  of  the  incompressible  equations  depends  on  the  ini¬ 
tial  conditions  through  the  parameters  ( SK/e)0  and  Re\  0.  The  compressible  case  depends 
additionally  on  the  initial  levels  of  the  thermodynamic  fluctuations,  dilatational  velocity  com¬ 
ponent  and  turbulent  Mach  number.  Apart  from  the  influence  of  initial  conditions,  another 
and  perhaps  more  important  feature  is  the  influence  of  local  compressibility  and  interactions 
between  the  velocity  and  thermodynamic  fields  on  the  statistics  at  a  given  time.  Equilibrium 
scalings  (if  present)  are  also  important  to  ascertain  so  as  to  improve  the  present  capabilities 
of  modeling  compressible  turbulence.  In  order  to  address  these  issues,  we  present  selected 
results  on  second-order  moments,  pdf’s  (probability  density  functions),  and  higher-order 
moments  of  both  the  velocity  and  thermodynamic  fields. 

3.1  Second-order  moments 

Evidence  from  physical  and  numerical  experiments  (  e.g.,  Tavoularis  and  Corrsin  (1981)  and 
Rogallo  (1981)  )  indicate  that  for  incompressible  homogeneous  shear,  both  the  turbulent 
kinetic  energy  and  the  turbulent  dissipation  rate  increase  exponentially.  Though  our  sim¬ 
ulations  are  in  accord  with  this  picture  of  exponential  growth,  we  note  that  Bernard  and 
Speziale  (1990)  have  recently  proposed  an  alternative  picture  where,  after  a  substantially 
long  period  of  exponential  growth,  the  turbulent  energy  eventually  becomes  bounded  for 
large  St  through  a  saturation  induced  by  vortex  stretching.  In  order  to  explore  the  effect 
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of  initial  levels  of  compressibility  on  the  evolution  of  the  turbulent  statistics,  we  have  per¬ 
formed  simulations  based  on  two  distinct  types  of  initial  conditions.  The  simulations  of  the 
first  type,  whose  results  are  shown  in  Figs.  3a-6a,  start  with  incompressible  data,  that  is, 
Prm»,o  =  Xo  =  0,  but  hav„*  different  initial  Mach  numbers  Mt,o-  The  second  type  of  simula¬ 
tions,  whose  results  are  shown  in  Fig.  3b-6b,  start  with  different  levels  of  initial  rms  density 
fluctuations  rp0  —  (Prm*/p)o  and  compressible  fraction  of  kinetic  energy  Xo,  but  have  the 
same  Aftto  =  0.3.  We  note  that  all  the  cases  of  Figs.  3-6  have  Rex, o  =  24  and  (SK / e)0  =  7.2. 
The  DNS  results  of  Fig.  3a  show  that  the  level  of  the  Favre- averaged  kinetic  energy  K  at 
a  given  time  decreases  with  increasing  Mt> o  in  the  simulations  with  the  first  type  of  initial 
conditions,  while  Fig.  3b  shows  that  the  level  of  K  also  decreases  with  increasing  Xo  and  rPy 0 
for  the  simulations  with  the  second  type  of  initial  conditions.  Thus  an  increase  in  compress¬ 
ibility  level,  either  due  to  increased  Mach  number  or  increased  dilatational  fraction  of  the 
velocity  field  decreases  the  growth  of  turbulent  kinetic  energy  in  the  case  of  homogeneous 
shear  flow.  Fig.  4  shows  that  the  development  of  the  turbulent  dissipation  e  has  two  phases. 
In  the  first  phase  (St  <  4)  the  higher  compressibility  cases  have  higher  values  of  e  probably 
due  to  the  buildup  of  the  compressible  dissipation  ec.  Later,  in  the  second  phase,  the  trend 
reverses  and  the  compressible  cases  have  lower  e  relative  to  the  incompressible  case.  Thus, 
compressibility  results  in  decreased  growth  of  both  e  and  K . 

In  order  to  explain  the  phenomenon  of  reduced  growth  rate  of  kinetic  energy,  we  consider 
the  equation  governing  the  kinetic  energy  of  turbulence  in  homogeneous  shear  which  is 

Jj.(pK)=  -Spu'{v%  -pe  +  p'd'  (8) 

where  pe  is  the  turbulent  dissipation  rate  and  p'd'  is  the  pressure-dilatation.  The  overbar 
over  a  variable  denotes  a  conventional  Reynolds  average,  while  the  overtilde  denotes  a  Favre 
average.  A  single  superscript  '  represents  fluctuations  with  respect  to  the  Reynolds  average, 
while  a  double  superscript  "  signifies  fluctuations  with  respect  to  the  Favre  average.  It  was 
shown  in  Sarkar  et  al.  (1989)  that  the  effect  of  compressibility  on  the  dissipation  rate  can 


7 


be  advantageously  studied  by  using  the  decomposition 

e  =  e3  +  ec  (9) 

where  the  solenoidal  dissipation  rate  t,  —  uuj[u[  and  the  compressible  dissipation  rate  ec  = 
(4/3 )ud'2 .  Here  u>[  denotes  the  fluctuating  vorticity  and  d1  denotes  the  fluctuating  divergence 
of  velocity.  We  note  that  the  correlations  involving  the  fluctuating  viscosity  have  been 
neglected  on  the  rhs  of  Eq.  (9).  Substituting  Eq.(9)  into  Eq.(8)  gives 

^-(pK)  =  -Spu "u'i  -  pe,  -  pec  +  p'd'  (10) 

at 

and  the  last  two  terms  represent  the  explicit  influence  of  the  non-solenoidal  nature  of  the 
fluctuating  velocity  field  in  the  kinetic  energy  budget.  The  DNS  time  histories  for  ec  and  p'd' 
are  presented  in  Figs.  5-6  respectively.  The  compressible  dissipation  ec  is  always  positive  (as 
expected)  and  when  starting  from  a  zero  initial  value  as  in  Fig.  5a  shows  a  monotonic  increase 
with  time.  The  pressure-dilatation  p'd'  is  highly  oscillatory  as  seen  in  Fig.  6.  Though  p'd' 
can  be  both  positive  or  negative  it  tends  to  be  predominantly  negative.  The  tendency  of 
p'd'  to  be  negative  in  homogeneous  shear  is  clearly  shown  in  Fig.  7  where  the  oscillations  are 
smoothed  out  by  plotting  values  of  p'd'  obtained  by  averaging  over  successive  time  intervals. 
Interestingly  enough,  on  comparing  Figs.  5a  and  7a  it  seems  that  the  time  averaged  values 
of  -p'd'  and  ec  are  of  the  same  order  and  exhibit  similar  trends  of  increase  with  St  and  Mti0- 
These  DNS  results  along  with  the  kinetic  energy  equation  (10)  show  that  both  ec  and  p'd' 
have  a  dissipative  effect  on  the  turbulent  kinetic  energy. 

The  predominantly  negative  values  of  p'd'  in  homogeneous  shear  is  in  contrast  to  the 
predominantly  positive  values  reported  in  Sarkar  et  al.  (1989)  for  the  case  of  decaying 
isotropic  turbulence.  To  understand  this  contrasting  behavior  of  p'd'  we  write  the  exact 
equation  for  the  pressure  variance  which  is  applicable  to  both  these  flows 

j-p'2  =  -lypp'd'  -  (27  -  1  )p'p'd'  -  2cp  +  <pp  (11) 

at 

where  <f>p  is  a  term  depending  on  mean  viscosity  and  mean  conductivity  which  is  negligible 
compared  to  ep.  Here  tp  -  (7  -  l)RKpT'}T'}  and  can  be  called  the  pressure  dissipation 
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term  by  virtue  of  being  a  sink  on  the  rhs  of  Eq.  (11).  Comparing  Eq.  (10)  and  Eq.  (11)  it 
is  clear  that  p'd'  acts  to  transfer  energy  between  the  kinetic  energy  of  turbulence  pK  and 
the  potential  energy  of  the  turbulence  p'2  j{2yp).  In  the  case  of  homogeneous  shear  the  rms 
values  of  the  velocity  and  pressure  increase  with  time.  Now,  the  third  term  on  the  rhs  of 
Eq.  (11)  is  always  negative  and  therefore  a  sink  for  the  pressure  variance.  For  small  p'/p  we 
can  neglect  the  second  term  on  the  rhs  of  Eq.  (11)  with  respect  to  the  first  term.  Therefore 
for  p'2  to  increase  with  time,  a  source  term  is  necessary  which  implies  that  p'd.1  be  negative 
in  accord  with  the  DNS  results.  We  note  that  a  similar  analysis  of  the  equations  for  p'2  and 
T'2  indicate  that  p'd'  and  T'd'  also  have  to  be  predominantly  negative  in  this  flow. 

In  decaying  isotropic  turbulence  p'2  has  to  decay  with  time;  however,  because  ep  is  suffi¬ 
cient  to  ensure  decay  of  p'2  the  sign  of  p'd' cannot  be  determined  from  Eq.  (11).  Alternatively, 
we  consider  the  exact  density  variance  equation  which  for  homogeneous  shear  flow  is 

-j-p'2  =  —2pp'd'  —  p'p'd'  (12) 

dt 

For  small  p' / p,  the  first  term  on  the  rhs  of  Eq.  (11)  dominates  the  second  term.  Therefore 
for  p'2  to  decrease  with  time  it  is  necessary  that  p'd'  be  predominantly  positive.  If  we 
assume  that  the  thermodynamic  fluctuations  are  approximately  isentropic  in  this  case  (DNS 
supports  this  assumption),  it  immediately  follows  that  p'd'  is  also  predominantly  positive 
in  decaying  isotropic  turbulence.  Thus,  the  role  of  p'd'  as  a  mechanism  for  energy  transfer 
between  the  kinetic  energy  and  potential  energy  dictates  its  differing  signs  in  homogeneous 
shear  and  unforced  isotropic  turbulence. 

Fig.  8a  shows  the  temporal  evolution  of  Mt  for  different  choices  of  the  initial  values  for 
Mt  and  Rex .  For  all  the  cases,  the  turbulent  Mach  number  shows  a  monotonic  increase 
with  time  for  the  range  of  St  spanned  by  the  simulations.  Two  results  of  the  analysis  of 
Sarkar  et  al.  (1989)  which  were  validated  previously  in  isotropic  turbulence  are  now  checked 
in  homogeneous  shear  flow.  These  results  are,  first,  the  compressible  dissipation  obeys  the 
scaling  ec/e,  —  aM2,  and  second,  the  partition  factor  F  ~  1  (see  Eq.  (7)  for  the  definition  of 
F )  which  implies  a  tendency  toward  equipartition  between  the  kinetic  and  potential  energies 
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of  the  fluctuating  compressible  mode.  Fig.  9a  shows  that  for  a  variety  of  DNS  cases  the 
ratio  ec/(e,M?)  approaches  an  equilibrium  value  of  approximately  0.5  for  large  enough  St. 
The  same  ratio  is  plotted  as  a  function  of  local  Mt  in  rig.  9b.  Thus  Fig.  9  supports  our 
previously  derived  result  on  the  dependence  of  the  compressible  dissipation  on  turbulent 
Mach  number.  Zeman  (1989)  had  proposed  the  scaling  eje,  —  1  —  exp{[(A/t  —  0.1)/0.6j2} 
which  is  clearly  in  contradiction  with  the  DNS  results  of  Fig.  9.  The  partition  factor  F  is 
seen  in  Fig.  10a  to  approach  and  oscillate  around  an  equilibrium  value  of  approximately  0.95 
indicating  that  equipartition  in  the  energies  associated  with  the  compressible  mode  holds  in 
the  case  of  homogeneous  shear  too.  In  Fig.  10b  we  show  the  behavior  of  Ft  =  7 ,!M?x/p'2 
which  is  defined  using  the  total  pressure  variance  rather  than  the  compressible  pressure 
variance  appearing  in  the  definition  of  F.  The  parameter  Ft  also  appears  to  reach  values 
independent  of  initial  conditions,  however  this  trend  is  less  pronounced  in  Ft  relative  to  F. 

In  Fig.  11  we  shew  the  time  evolution  of  rms  density  and  temperature  fluctuations.  It 
appears  that  both  the  rms  normalized  density  and  temperature  fluctuations  eventually  scale 
as  A/2.  The  rms  pressure  also  scales  as  A/2.  Erlebacher  el  al.  (1990)  decomposed  the  pressure 
p  into  an  incompressible  component  p 1  and  a  compressible  component  pc .  The  beho.  ior  of 
P?mt  an^  pimt  shown  in  Fig.  12.  The  instantaneous  incompressible  pressure  p1'  satisfies 
the  Poisson  equation  Eq.  (6).  The  compressible  pressure  pc'  =  p'  —  p1'  is  the  remaining 
part  of  the  pressure  which  arises  from  both  density  variations  and  non- zero  divergence  of  the 
velocity.  Fig.  12  indicates  that  the  rms  of  both  components  of  the  normalized  pressure  vary 
like  A/2.  A  simple  order  of  magnitude  analysis  of  Eq.  (6)  along  with  recognition  of  the  fact 
that  SKje  —  0(1)  shows  that  plnaBjp  —  Thus  it  seems  that  in  homogeneous  shear 

the  conventionally  accepted  Mach  number  variation  holds  for  pIrms-.  p^,  and  pnn9. 

The  evolution  of  pc  in  the  case  where  the  initial  data  is  incompressible  can  be  understood 
by  considering  the  following  equation  for  compressible  pressure  (see  Erlebacher  et  al.  (1990) 
for  details)  where  heat  conduction  and  viscous  effects  hav  been  neglected 


C f  I '  c f 

o,p  +  u,  p  '  , 


C  C' 


-1  V 


,  Cl' 
1-  ul  p 


f  7(  1  1  pC') 


c 1 


u. 


(dtP! 


V  /' 
+  Ux  P 


(13) 
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Since  pc'  and  uf '  are  initially  zero  the  growth  of  pc'  is  due  to  the  forcing  by  the  incompressible 
pressure  on  the  rhs  of  Eq.  (13).  This  forcing  causes  p^/p  to  become  0(prrms/p)  —  O(Mf)  on 
the  convective  time  scale  which  explains  the  trend  seen  in  the  three  lower  curves  of  Fig.  12b. 
The  top  curve  of  Fig.  12b,  which  corresponds  to  a  case  with  initially  non-zero  p^  and  u^, 
also  shows  the  same  trend  indicating  that  the  forcing  by  the  incompressible  velocity  field 
(through  p1  )  has  a  substantial  influence  in  determining  the  level  of  p^. 


3.2  Probability  density  functions  and  higher  order  moments 

The  probability  density  function  (pdf)  gives  a  comprehensive  description  of  a  fluctuating 
variable  because  it  formally  enables  the  calculation  of  all  the  statistical  moments.  Second- 
order  moments  (discussed  in  the  previous  section)  are  probably  the  most  important  statistical 
indicators  of  the  turbulence  because  they  characterize  intensities  and  cross-correlations  of  the 
fluctuating  variables.  Higher-order  moments  such  as  the  skewness  and  flatness  give  additional 
information  on  the  statistical  distribution  of  the  variable,  for  example,  a  positive  skewness 
indicates  a  pdf  which  is  skewed  towards  positive  fluctt  .tions  from  the  mean,  while  a  large 
flatness  indicates  a  large  propensity  for  the  occurrence  of  extreme  values  of  the  variable.  The 
skewness  Sk(9 )  and  flatness  Fl{9)  of  a  variable  9  are  defined  by 


Sk(9)  = 


Fl(9)  = 


9* 

(P)3/2 

0* 

(Py 


A  quantity  with  a  Gaussian  pdf  has  Sk  =  0  and  FI  =  3.  In  this  section,  we  present  results  on 
pdf,  skewness  and  flatness  of  some  variables  for  a  DNS  case  with  significant  compressibility 
levels.  The  initial  nondimensional  parameters  for  this  case  are 


Re\  =  24  ,  SK/e  =  7.2  ,  Mt  =  0.4  ,  *=0  ,  Prm/p  =  0  ,  p^/p=0 


The  pdf  /„  of  the  streamwise  fluctuating  velocity  in  Fig.  13a  appears  to  be  symmetric 
around  the  origin.  The  skewness  and  flatness  of  u  have  insignificant  deviations  from  Gaussian 


values.  The  experiments  of  Tavoularis  and  Corrsin  (1981)  in  homogeneous  shear  flow  and 
Frenkiel  and  Klebanoff  (1973)  in  the  log-region  of  the  flat  plate  turbulent  boundary  layer 
also  indicate  that  u  is  approximately  Gaussian.  The  pdf’s  of  the  other  components  of 
the  velocity  are  very  similar  to  those  of  the  streamwise  velocity  component.  However  the 
transverse  velocity  v  develops  a  somewhat  higher  flatness  ( Fl(v )  =  3.3  at  St  =  8  )  than 
the  other  velocity  components.  Fig.  13b  shows  the  pdf  f m  of  the  Mach  number  M  =  q/c 
where  q  is  the  instantaneous  flow  speed  and  c  is  the  instantaneous  speed  of  sound.  M  has 
maximum  values  of  1.7  and  3.2  at  St  =  1  and  St  =  8  respectively,  indicating  that  the  pdf  of 
M  has  exceptionally  long  tails. 

We  now  investigate  quantities  related  to  the  velocity  derivatives.  Figs.  14a-b  show  the 
pdfs  of  the  spanwise  vorticity  component  uz  and  the  dilatation  d.  Both  pdfs  are  skewed 
towards  negative  values.  Higher-order  moments  of  the  vorticity  components  and  dilatation 
are  shown  in  Fig.  15.  The  skewness  of  the  dilatation  (see  Fig.  15a)  reaches  values  of  about 
—0.5,  signifying  that  local  compressions  are  more  likely  than  local  expansions  in  the  flow. 
Because  of  the  invariance  of  the  homogeneous  shear  problem  to  a  reflection  of  the  spanwise 
(z  direction)  axis,  the  u)x  and  u>v  components  have  zero  skewness.  The  dilatational  field 
has  larger  flatness  than  the  vorticity  field,  indicating  that  it  is  more  intermittent  than  the 
vorticity  field. 

In  this  DNS  case,  the  maximum  magnitude  of  dilatation  in  the  compression  regions 
attains  values  which  are  about  an  order  of  magnitude  larger  than  the  rms  dilatation.  Lee, 
Lele  and  Moin  (1990)  have  also  observed  such  high  values  of  the  maximum  dilatation  in 
isotropic  turbulence,  and  have  furthermore  associated  these  high  values  with  dynamically 
significant  eddy  shocklets.  However,  in  order  to  assess  the  importance  of  these  high  dilatation 
regions  it  is  necessary  not  only  to  identify  the  magnitude  of  the  maximum  dilatation  but 
also  to  calculate  the  contribution  of  these  regions  to  the  rms  dilatation.  We  find  that,  at 
St  =  8,  the  compressive  regions  with  a  magnitude  of  dilatation  larger  than  3  times  the  rms 
dilatation  contribute  only  about  2%  to  the  variance  of  dilatation.  This  suggests  that  shocklet 
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dissipation  is  not  the  dominant  contributory  mechanism  to  the  compressible  dissipation  for 
Mt  <  0.5. 

From  Fig.  16a  it  appears  that  the  imposition  of  shear  does  not  alter  the  skewness  and 
flatness  of  the  Mach  number  M  appreciably  from  the  isotropic  turbulence  values  at  St  =  0. 
The  fact  that  M  is  bounded  from  below  and  not  from  above  may  be  responsible  for  the 
departure  of  skewness  and  flatness  from  Gaussian  values.  The  instantaneous  values  of  M 
can  become  rather  large  (  the  ratio  M^x/Mt  =  5.8  at  St  —  8),  but  me  flatness  level  of  4 
at  St  =  8  indicates  that  these  large  extrema  have  a  relatively  low  probability  of  occurrence. 
The  long  tail  in  the  pdf  of  M  (see  Fig.  13b  and  the  previous  discussion  of  the  figure)  does  not 
seem  to  contribute  significantly  at  the  level  of  third  and  fourth  moments.  Fig.  16b  shows  that 
the  velocity  field  can  be  highly  supersonic  =  3.3  at  St  =  8).  Though  the  percentage 

of  the  volume  occupied  by  supersonic  fluctuations,  denoted  by  <^,uper  in  Fig.  16b,  increases 
with  time,  its  maximum  value  is  a  modest  4.9%  at  St  =  8. 

4  Characteristics  of  the  instantaneous  flow 

The  previous  sections  have  focused  on  the  influence  of  compressibility  on  the  statistical 
attributes  of  homogeneous  shear  turbulence.  We  now  present  some  preliminary  results  on 
the  spatial  characteristics  of  the  instantaneous  turbulence  from  ongoing  flow  visualization. 
Figs.  17-20  are  results  from  the  flow  field  at  St  =  9  for  a  case  which  has  the  following  initial 
nondimensional  parameters 

Rex  =  21  ,  SK/e  =  5.8  ,  Mt  =  0.3  ,  x  =  0.09  ,  p^/p  =  0.09  ,  p^/p  =  0.13 

By  St  =  9  the  initially  isotropic  turbulence  develops  into  a  realistic  representation  of  ho¬ 
mogeneous  shear  turbulence.  The  turbulence  at  St  =  9  has  Mt  =  0.42,  prmm/p  =  0.12  and 
X  =  0.06,  which  implies  a  moderate  level  of  compressibility,  and  has  a  Rex  =  32.  One  of 
the  important  characteristics  of  the  velocity  field  is  the  extent  and  topology  of  supersonic 
regions.  Fig.  17  shows  a  3-D  contour  plot  of  regions  with  fluctuating  Mach  number  M  >  1. 
These  regions  of  supersonic  fluctuations  are  small  discrete  regions  scattered  through  the  flow 
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domain,  and  are  elongated  in  the  streamwise  direction.  The  volume  occupied  by  the  regions 
with  supersonic  fluctuations  is  small.  Fig.  18  shows  the  dilatation  V.u'  in  a  plane  parallel 
to  the  (xi,x2)  plane.  The  regions  in  black  are  expansion  regions  with  positive  dilatation, 
while  the  regions  in  grey  and  white  are  compression  regions  with  negative  dilatation.  Both 
the  compression  regions  and  the  expansion  regions  are  very  elongated  in  the  streamwise  xx 
direction.  The  white  portions  in  the  dilatation  plot  of  Fig.  18  correspond  to  compression 
regions  where  the  magnitude  of  dilatation  is  larger  than  the  maximum  positive  dilatation. 
The  magnitude  of  dilatation  in  these  white  regions  are  up  to  50%  larger  than  the  maximum 
positive  dilatation  indicating  that  compressions  are  stronger  than  the  expansions.  As  Fig.  18 
shows,  the  white  regions  of  strong  compression  are  small  and  discrete  elements  which  occupy 
a  small  portion  of  the  plane.  Figs.  19a  and  19b  are  contour  plots  of  the  magnitude  of  u\'  and 
respectively  in  a  plane  parallel  to  the  (xx,X2)  plane.  Regions  with  large  gradients  seem 
to  be  more  prevalent  in  the  uf  field  relative  to  the  u\'  field.  Figs.  20a  and  20b  show  the 
incompressible  pressure  p1'  and  the  compressible  pressure  pc'  in  the  same  plane  considered 
in  Fig.  19.  It  is  clear  that  the  p1'  field  has  significant  differences  with  respect  to  the  pc>  field. 
The  compressible  pressure  has  higher  gradients  and  more  pronounced  small  scale  features 
than  the  incompressible  pressure.  Elongated  and  inclined  structures  are  also  visible  in  the 
pc'  field. 

5  Conclusions 

We  have  performed  direct  numerical  simulations  of  homogeneous  shear  flow  using  a  spectral 
collocation  technique  and  analyzed  the  data  base  in  order  to  quantify  compressibility  effects 
on  the  turbulence.  An  increase  in  the  compressibility  level  (i.e.  turbulent  Mach  number, 
density  fluctuations  or  rms  dilatational  velocity  fluctuations)  leads  to  a  decrease  in  the 
growth  of  turbulent  kinetic  energy.  In  the  compressible  case,  the  kinetic  energy  equation 
has  two  additional  terms,  the  compressible  dissipation  and  the  pressure-dilatation.  Both 
these  terms  have  a  dissipative  contribution  in  shear  flow,  leading  to  the  reduced  growth  of 
the  turbulent  kinetic  energy.  The  0(Mf2)  variation  of  the  compressible  dissipation  which 
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was  obtained  previously  by  an  asymptotic  theory  and  then  verified  in  the  DNS  of  isotropic 
turbulence,  has  now  been  substantiated  by  the  DNS  of  homogeneous  shear  turbulence.  The 
pressure-dilatation  is  a  more  difficult  term  to  characterize,  because  it  can  be  either  positive 
or  negative  and  also  because  gradient  transport  arguments  are  clearly  inapplicable.  Our  DNS 
show  that  the  pressure-dilatation  is  predominantly  negative  in  the  homogeneous  shear  case 
in  contrast  to  its  predominantly  positive  character  in  decaying  isotropic  turbulence.  Analysis 
shows  that,  in  the  transport  equation  for  the  pressure  variance,  the  pressure-dilatation  has 
to  be  a  source  in  the  case  of  homogeneous  shear  (where  the  turbulence  grows)  and  a  sink 
in  the  case  of  decaying  isotropic  turbulence,  thus  accounting  for  the  opposite  signs  of  the 
pressure-dilatation  in  the  two  flows.  In  the  core  region  of  inhomogeneous  aerodynamic  flows, 
perhaps  the  pressure-dilatation  is  predominantly  negative  so  as  to  transfer  energy  from  the 
velocity  field  to  the  pressure  field. 

It  is  well  known  that,  for  low-speed  flows,  the  rms  pressure  normalized  by  the  mean 
pressure  varies  as  0(Mt2).  Our  DNS  results  support  this  scaling  of  the  rms  of  pressure, 
density  and  temperature  fluctuations  for  compressible  homogeneous  shear  flow.  These  DNS 
results  and  our  analysis  suggest  that  such  a  scaling  may  be  valid  for  those  inhomogeneous 
flows  where  the  thermodynamic  fluctuations  arise  due  to  high-speed  effects  and  not  due 
to  other  factors  such  as  external  heating  or  cooling,  and  mixing  of  fluids  with  different 
densities  or  temperatures.  We  have  studied  the  probability  density  functions  and  higher- 
order  moments  of  the  dilatation  and  the  vorticity.  The  dilatation  field  has  appreciable 
negative  skewness  indicating  a  preponderance  of  compressions  over  expansions,  and  has  a 
larger  flatness  than  the  vorticity  indicating  that  it  is  more  intermittent  than  the  vorticity. 
The  intermittent  regions  of  negative  dilatation  were  found  to  have  only  a  small  contribution 
to  the  compressible  dissipation  for  Mt  <  0.5. 
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Figure  lib  Behavior  of  scaled  rms  temperature. 


Figure  12a  Behavior  of  scaled  rms  incompressible  pressure. 


Figure  12b  Behavior  of  scaled  rms  compressible  pressure. 
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Figure  16a  Skewness  and  flatness  of  Mach  number. 


Figure  16b  Maximum  Mach  number  and  volume  fraction  of  supersonic  fluctuations 


Figure  19a  Contours  of  the  magnitude  of  u[  in  a  (11,12)  cut  of  the  flow  domain. 


Figure  19b  Contours  of  the  magnitude  of  uf 
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Figure  20a  Contours  of  p1  in  the  plane 
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